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SYNOISIS 

AH AHAIirTICAI STUDY 01’ THE THBEMAL BEIECTS OP HEAT HEJECTIOH 
PROM POV®R PIAHTS TO COOI1III& POISIB 

A Thesis Submitted 

In Partial Pulfilnent of the Requirements 

Por the Degree of 

MASTER OP TBCHHOID&y 
i3y 

SHEOIEE TItAEY 
to the 

Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur 
July 1978 

In the present work, the steady state and transient, verti- 
cal temperature distribution in a cooling pond have been deteimined 
analytically. The effect of various governing parameters, such as 
ambient temperature, cloudiness, relative humidity, wind velocity, 
solar flux, eddy diffusivity and power plant capacity on the tempera- 
tuxes are studied in detail for a one dimensional model. The govern- 
ing equations are solved numerically by finite difference meihod « A 
two-dimensional model, for the limiting cases of inviscid and creep- 
ing flow, neglecting variations in the lateral direction is also studied 
to determine the recixoulation in the cooling pond. The study deter- 
mines the effect on the intake temperature due to the heat rejected 
by a power plant into a water body* 



CHAPTER 1 


IKE'RODUOTIOl 


1.1 The Problem : : 

Electric power has become a virtual necessity for the healthj 
ccmfort and economic well being of the people. At present ja major 
part of the electric energy produced is generated by steam driven 
thermal power plants. Because of the potential impact on the quality 
of the air and water environment, due to heat rejection, the problem 
of setting up a£ large, new thermal power plants is beccming increas- 
ingly critical. One of the most crucial considerations is that of the 
disposal ctf the waste heat, inherent in the generating process keeping 
the efficiency of the energy generation system high and the effect on 
the environment low. 

A characteristic of the operation of a steam driven thermal 
power plant is that a large flow of water is required in the condenseis 
to condense the exhaust steam, from the turbine, as water. This keeps 
Uie temperature of condensatian lowj so as to maximize the energy 
ccnversion. The condensate is recirculated back to the boiler. In 
flowing through the condensers, the cooling water is genesrally heated 
by 6°c to 9**o, depending upon the design of the plant. The amount 
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of waste heat discharged to the condenser is directly related to the 
efficiency of energy genei^aticn of the plant. An indication of the 
magnitude of the waste heat is the fact that even at the most effi- 
cient plants now in operatioa, the heat rejected is generally muda 
greater than the heat equivalent of the electric energy generated. 

At atypical efficiency of 33^ '^-twice the amount of energy generated 
is rejected as waste heat, 

The heat added to -the water as it flows through the conden- 
sers must eventually he dissipated to the atmosphere. Many power 
plants use a cooling system for the dissipation of the waste heat, 
in which the water is tsfcen frcsn a river, lake, reservoir or the sea, 
passed throu^ the condenser and then returned at a higher 'tempeihattire 
to the water body. The added heat is suhsequentjy dissipated by means 
of increased evaporation, back radiation and conduction to the atmos- 
phere. The only consumptim of water is that resulting from an 
increased evaporation in the water body due to the addition of 
heat. 

At inland locations inhere natural streams or lakes of ade- 
quate size are not available and if suitable sites can be found, 
cooling ponds may be constructed to provide the cooling water needs'* 
The structural and flow requirements for cooling pohds are essentially 
the same as that for natural water bodies. However, it does entail 
an additional expenditure in constwction, \diich is of particular 
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significance for larger power plants# A natural attempt would j 
therefore |be to tap the natural resources such as lakes f rivers and 
the sea# 

Thn cooling water is circulated betv\reen the condensers and 
the water body# This generates a flow field in the \7ater body "vdiich 
results in a recirculation of the warm water, discharged at the outlet, 
to the cold water intake# This recirculation, which depends on the 
position of the intake and the outlet, the amount of heat rejected, 
flow rates, ambient conditions and the dimensions and characteristics 
of the water body, raises the temperature at the intaie and thus 
lowers the efficiency, of the pov/er plant, which is linked, from ther- 
modynamic considerations to the temperature at which heat is rejected, 
being lower for a higher temperature* 

Most of the water bodies , under Indian conditions are gene- 
rally stratified over a large part of the year# This implies that 
a heated upper layer lies over a colder layer# This is a stable situa- 
tion since colder layers are heavier than heated layers# Unless forced 
by restrictions on the outfall temperature frcm the power plant, the 
heated discharge is usually floated over the surface of the stratified 
water body, in order to obtain a high surface temperature, ’vdiich gives 
rise to a high rate of heat loss at the surface and minimum effect on 
the water body# At the same time, the intake is located below the 
interface between the two layers, termed as Ihe thermocline, so as 
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to employ the lowest available temperature oC cooling water, for 
increased efficiency of energy generation. A general sketch of the 


problem of recirculation is shown in Pig. 1. 

1 *2 Review of Previous Work ; 

Initially, this area was largely considered due to the en~ 
vircmental effects arising from heat re ejection. Recently, the 
need of a fluid flow and heat transfer analysis has been realised, 
in an attempt to improve the power plant efficiency. A comprehensive 
review of the literature with extensive references is given by Jaluria, 
Variyar and Mehta (1976). 

Raphael (1962) developed a procedure for predicting the tem- 
perature of various water bodies, such as shallow lakes, flowing 
streams and detention reservoirs, from weather records, inflow and 
outflow characteristics, the surface area and the volume of the water 
body. He assumed that the thermocline is absent and the water is so 
stirred by wind or internal currents that the temperature is unifona 
throughout. Relay and Seaders (1966) developed a mathematical model 
for predicting temperatures in rivers and reservoirs. They computed 
the net mcntiily energy exchange from the resejcvoir and Ihen distributed 
it vertically so that tile resulting temperature gradients approximated 
the -standard for each month. Rake and Horleman (1969) developed a 
theoretical model for the time dependent vertical tempera-ture distri— 
buticG in a deep lake, during -Khe yearly cycle. They took a heat-fl«x 
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balance at the water surface, vhich accounts for back radiation and 
evaporative loss as a boundary ccnditicn-. 

famai, ii?iegel and Tomberg (1969) realised that the main prob- 
lem, in the design of cooling water systems for themal power plants, 
is the determination of the nature of mixin.'^ of the warm water dis- 
charge, so as to find methods to minimize tb.e recirculation of the 

i 

waste heat. With this in mind, they performed several experiments, 
on the mixing of buoyant flows, discharged ) crizontally at the sur- 
face of a body of water, in the laboratory. They also collected data, 
from a number of sources , on the cooling waser capacities o£ thermal 
power plants, along wilh the flow characterj.stios , and compared them 
with the experimental values. Brown (1970) has discussed the various 
methods of waste heat disposal and has emphasised the need to develop 
improved and less costly cooling devices •F'^ef an (1970, 1972) showed 
that the flow and the cooling of a warm water surface layer can bo 
represented by or physical model, if the inf? ow and outflow of water 
and the atmospheric conditions, it is exposed to, are ccntrollable . 

He simulated the phencmena of mixing at liie outlet, stratification 
and -taie surface heat transfer through the use of physical models* 

He also described an experimental design to simulate a heated water 
discharge from a cbannel, into a deep lake for lateral spread* 

Woare and Jaluria (l97l, 1972) studied the thermal effects 
of power plants on lakes, in terms of its temperature cycle* They 
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developed a oiae-diiiiensional model based an the assumption that a 
Richardson number is constant at the base of any stratified layer'. 

This gave a constant heat flux across the theimocline. The model is 
then perturbed in terms of this heat fluxj as well as the thermal 
diffusivity to give the power plant impact on a given lake. They 
coaputed the change in the summer maximdn, winter minimum, stratifica- 
tion and overturn, or destratification, temperatures and time of occu- 
rrence .G-rubert and Abbott (1972) gave the mathanatical formulatiaa of 
nearly horizontal stratified flows, in terms of partial differential 
equations and their characteristics* They have also shown that the 
characteristic directions in a stratified fluid devide into pairs, 
each pair being associated with a fluid layer. Bindley and Mner 
( 1972 ) observed that the temperature distributicn in a cooling water 
body, being employed for a.' power plant, is dependent upon the manner 
in which the heated effluent mixes with the receiving water body and 
upon the rate of heat exchange between the water surface and the atmos- 
phere. They found that, far from the receiving end, where nilxing is 
nearly cmplete, the dominant cooling mechaniaa is heat exchaaoge wiUi 
the atmosphere. The rate of this cooling is directly proporticaaal to 
the swrf ace heat exchange co-efficient , which is in turn, a function of 
several environmental variables, particularly the wind veloci'ty'i 

Based upon 1iie assumption of horizontal isotherms, at all 
times, Huber, Harleman and Ryan (1972) developed a natbematicsl model, 
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for detennining the vertical temperature distributicaa in stratified 
reservoirs, udaich includes the effects of heat sources and sinks at 
the boundaries, internal absorpticn of solar radiation and heat trans- 
port b^ convection and diffosicn. Miyazaki (1974-) studied the heated 
two-dimensional jet discharge at the water surface. He carried out 
an analysis at vaidous Richardson numbers and calculated the velocity 
and temperature distributions. Snider and Yiskanta (1975) carried out 
an analysis for the time dependent thermal stratification in the surface 
layers of stagnant vsater due to solar radiation. Ihey used finite 
difference methods to obtain the transient temperature distribution, 
by solving the one-dimensional energy equation, for combined conduction 
and radiation energy transfer mechanisms. 

Jaluria, Variyar and Mehta (1976, 1977) developed a mathemati- 
cal model to predict the temperature and velocity profiles in a body 
cf water, due to the basic mechanisms in the presence of winds, cloudi- 
ness, back radiation etc. Ihey ccmputed the equilibrium average sur- 
face temperature, which is defined as the temperature attained by the 
surface, if the ambient conditiais are held constant at si»ciflc values, 
and studied the heat transfer medianism at the surface in detail. They 
also developed a more generalized model, considering convection effects 
caused by the outflow and intake cf cooling water and the heat addition 
from the power plant. They studied certain specific, s impli fied , caie 
and two-dimensional models in detail. The work concentrated largely 
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^)een 

ci£i steady state models and not much work has^done on the transient 
behaviour, which forms an Important consideraticn in the present 
study, as outlined below. 

1.3 Present Work : 

In the present work, a one-dimensicaial mathematical model, 
in the foim of fin±te difference equaticns, has been developed, to 
detenaine the vertical temperature distribution in the cooling pond# 

Both the steady state and transient, or time dependent, cases are 
studied in detail, for a natural lake, as well as for a lake, or 
cooling pond, with heat rejectim from power plants* 

The net surface heat exchange is computed by using a method 
similar to that outlined by Jaluria, Variyar and Mehta (1976) mainly 
for Indian conditions, with respect to the ambient conditions, such as 
solar flux, wind speed , ambient temperature, etc# The equilibrivaa 
average surface temperature is determined, which is employed in the 
computation of the vertical temperature distribution# The effects 
of various governing parameters such as ambient temperature, cloudiness, 
relative humidity, back radiation, soler flux, heat discharge from 
power plants, flow rates etc# on liiese results are studied in detail, 
to provide the first step in the design of such systems# 

The work considers 14ie transient temperature variation in 
detail, developing suitable models for stratification and destrati- 
fication of the water body# The effect of the power plant on the 
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intake teciperature rise is detemined and the dependence cxf this effect 
on the various parameters is outlined* It is found that the water body 
is stratified for most part of the year* The power plant tends to 
raise the temperature level and chaiige the stratification cycle of 
the water body. The effect of turbulence, is considered on the basis 
of the eddy diffiisivily and this is taken as a parameter, whose effect 
is also studied* The work is of considerable importance in the design 
of cooling water systems, with respect to flow rates, heat input, 
locations of intake and outfall, etc* 

A two-dimensional model, neglecting variations in the lateral 
direction, is also developed to detemine the recirculation in the 
cooling pond « The model is studied for the extreme cases of inviscid 
and creeping flow* The resulting velocity field is deteimined and 
the general nature of the flow ^ discussed* 

All the ccoputations were canned out cn IBM 7044 Computer, 
at IIT/Kanpur* The various meteorological data needed for the study 
were obtained largely from Indian Ifeteorological Department and partly 
from the analysis given by Raphael (1962)* 



CHAPTER 2 


THEOBETICAI. MD I EIS AlilD METHOD OP SOIUTIOR 

The behaviour of a water body^ which acts as a cooling water 
source for a power plant, is deteimnedjto a large extent, by the 
nature of its interaction with the environment, to which it ultimately 
dissipates the heat added by the power plant . Hence , in developiaig 
an analytical model for a ooolirg pond, the mathKnatical representa- 
tion of the various medianisms underlying its interaction wi13a Ihe 
envircaament is of considerable importance* The basic features of an 
analytical model for the lafce may now be outlined as follows: 

a) A finite pond or lake. 

b) Ho energy transfer across its sides and bottcm i.e., energy 
transfer "occurs only at its surface and at the intake and outfall 
for the cooling water of 1he power plant* This assumption is 
adequately supported by various investigaticns. 

c) An idealized configuration of vertical sides and flat horizaatal 
base would be considered initially* This could be eventually 
modified to consider, as close as possible, the exact topology 
of a given lake. 

d) The Outfall at the surface , Tfihile the intake at the bottom cxf the 
lake. This implies that the outfall enters the lake wilh the 
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smallest possible disturbance, This is applicable in most prac- 
tical applicaticsna . 

e) Temperature variatim in the natural lake is only in the vertical 
direction, the horizontal gradients being neglected, 

f) Ambient conditions are functions of time, though they may be 
teOcen as specific values over small time periods during the 
year. 

2,1 E nergy Exchange at the Surface : 

The energy exchange at the surface, is due to various heat 
transfer mechanisms outlined below; 

a) Radiation from the Sun an d from the Atmosphere ; 

At a point outside the atmosphere of the earth, radiant 
energy is received from the sun, on a surface normal to its rays, at 
a rate of approximately 1165 K.cal/hr/sq JU , . This value, called 

the ’^olar Oonstant”, fluctuates sli^tly during the course of the 
year because of sunspots and variation in the distance between the 
earth and the sun, A plane, normal to the rays of the sun at the 
surface of the earth receives considerably less energy, as compared 
to the solar constant because moisture, gases and solid particles in 
the air scatter and absorb much of the incident radiation, Althou^ 
much of the direct radiant energy is lost, some of this energy is 
recovered as dif:^e radiation. 
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Of "the 1:01:81 quantity of the direct, or short wave, radiation 
that reaches the surface of a lake, a pcKctiori is returned imcharged 
due to reflection, at the surface and scattering by bubbles and sus- 
pended particles immediately below it. However, for practical engin— ’ 
eering computations, solar radiation and reflected radiation can be 

combined in one functicn, (Q. - Q ) or s 

* i r 

Q (effective absorbed solar radiation) = Q. - Q' (2,1a) 
s 1 r 

If 9 c 9 is the average cloud cover, in tenths of sky covered, the net 

short-wave radiation is obtained from the correlation given by « 

Qg , ( 1 - 0.0071 c^) (Q^ - Q^) (2.1b) 

b) B ack Radiation ; 

Effective back radiation may be defined as the difference 
between , the lo33gwave radiation leaving a body of water, and 

longwave radiation from the atmosphere being '..absorbed by 
the body of water. Atmospheric radiation does not follow any simple 
law as it is a function of many variables, such as the distribution 
of moisture, temperature, oaone, carbon dioxide, etc. Thus, this has 
been found to be dependent on cloud hei^t, as hi^, middle and low 
clouds, and the amount of cloud, as scattered cloud, broken cloud and 
over cast. However, as most -v^eather observations give cloud amount 
in tenths of sky obscured, the variation of the atmospheric radiation 
factor has been given as a function of cloud amount and vapour 
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pressure by Eapliael (1962), The effective back radiatian, can 
be computed as f 

where, 0“is Stefan-Boltzman Constant, f is absolute water surface 

s 

temperature, is absolute ambient temperature and ^ is sky radia- 
tion factor, tabulated by Raphael (1962), 

c) Convectiaia ; 

Sensible heat is oonvected to or from a body of water when- 
ever a temperature difference exists between air and water. The basic 

equation for the canvectian heat transfer, Q , from the water surface 

c 

to the air is, 

= h (Tg - T^) (2,3) 

vhere h is the surface heat transfer coefficient. The value of 
h is largely dependent on the condition of 1he air in contact with 
the water and hence it is a functicn of the wind velocity, its tempe- 
rature, direction and turbulence level, etc. The following corre- 
lation, by Wada (1968), has been found to be most suitable for condi- 
tions relevant to India. 

h = 2,77 X lo“^ (0.48 + 0.272 v) K.cal/m^.sec. ^c (2.3a) 
■vdaere Y is the wind velocity along the water surface in meter/sec'. 

d) Evaporation : 


Yhen the vapour pressure of ambient air is less than the 
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saturated vapour pressure y at the water surface temperature j water 

evapcsrates into the air, removing heat frcaa the water mainly due to 

the energy reqtuired to evaporate the water and to a small extent, by 

the secondary heat transfer due to the sensible heat contained in the 

v/ater removed by evaporation. The correlation used to calculate the 

heat transfer by evaporation is obtained from the work of Bindley and 

Miner (1972). !Ehis is given as, 

Q = 2h Pe (T ) - e (T )} (2.4) 

e L s aJ 

where, e(0} ) is the partial vapour pressure at water surface tempera- 
s 

ture in mm of hg, ®(^g^) 1-® partial vapour pressure at ambient 

temperature in mm of hg and h is the heat transfer coefficient 

defined by Bqn. (2.5a). Also, 

6(1! ) = e (T ) X EJH. (2*4a) 

a s a 

where, e (t ) is the saturated vapour pressure at ambient temperature 
s a 

and E.H. is the relative humidity at ambient temperature. 

The above relationship can also be obtained by considerong 
14ie analogy between heat and mass transfer and the corresponding 
Prandtl Schmidt numbers. Therefore, the heat loss from the sur- 
face is given by , 

Q = Q +0, -Q„ + Q 

e br s c 

The surface equilibrium temperature T^, may often be repre- 
sented approximately as a sinusoidal variation over the year. This 
approKlmation has been used extensively by many investigators, see, 



15 


for instaaacej the represcntatiaa of Moore and Jaluria (1972) as* 


T = A - B cos 
e e e 


L 365 J 


( 2 '. 6 ) 


(t ) + (t ') 

where, . e'^max* e min. 
%= 2 


B = 


(T ) - (T ) . 

e max. e nnn. 


(2, 6a) 
(2.6b) 


and C is time in days at which the minimum surface equlibrium tempe— 

0 

rature occurs, where t is time in days. 

OJhe rate of net heat loss from the surface, Q, is given as* 

Q = K (T - T )■ (2*7) 

^ s e' 

i/daerei K is the surface heat exchange coefficient and the surface 

equilibrium temperature represented as a sinusoidal wave# In the 

above equation^ K and T are dependent on ambient conditions and 

e 


may be deteimined for a given lake by a detailed study of Q and 
in temffi of various mechanisms , as outlined above. 


2.2 One Dimens i goal Models* 

Ihysically, a one dimensional model represents the circums- 
tance when the lateral mixing in the lake, or the cooling pond, is so 
vigorous that only vertical exchange of heat and mass need be consi- 
dered. It has been found by Raphael (1962) that the eddy diflUsivily 
in the vertical direction, for stratified water bodies, is around one- 
hundredth the value in the horizontal direction. Therefore , horizontal 
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mixing is much more vigorous , mainly due to the huo 3 rancy forces 
resulting from the temperature dependent density differences. A uni- 
form temperature in the horizontal plane may be assumed over a finite 
area, of the water body, vfliere most of the flow occurs* A ccnvec- 
tive loss from the side to tae surrounding water medium, at the tem- 
perature for the natural lake, may be assumed. In actual practice, 
the recirculation would essentially occur over a portion of the water 
body and energy would be lost to the ne labouring water through a 
similar induced moticn, arising due to the entra in ment requirements 
of the recirculating flow. 

For a one dimensional model, since temperature variation 
exists only in the vertical direction, the boundary conditiona have 
to be specified at the surface and at the bottom, with the convec- 
tive loss, if any, to the neighbouring fluid, taken as an energy 
loss from each fluid element. The theoretical model is developed 
for the steady state, as well as for, the transient, as outlined 
below: 

2.2.1 Steady State : 

A side loss model with energy input from a power plant is 
considered. ITeglecting density variations , the continuity and energy 
equations yield the following, lixe momentum equation being only the 
hydrostatic pressure variation equaticn: 


dw 

dz 


0 


( 2 . 8 ) 
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r d T dC 

C „ ” w — 


vd-th the conditicxQS ! 


h (T - T ) 
w e 

C 1 
P 


= 0 


OP = I and w ss w„ at z « 0 
s 0 


T = T, and w = w. at z » H 

D 0 


(2.9) 

(2,9a) 

(2.9b) 


W 


+ a t) - + C 


dl ^ 
H dz 


= 0. at su3?face 


... (2.9c) 

where 5 w is- the vertical velocity in the cooling pend, ^ g is the 
eddy diffusivity, h^ the side convective heat transfer coefficient, 
0P.j^ the bottom element temperature of the cooling water body, /I T the 
difference between outlet and inlet temperature of condenser cooling 
water, I the side length of the cooling water body over vhich the 
temperature is taken as unifemm and A , the top surface area of the 
cooling pond. Here A T is. taken as a known quantity, h^ is varied 
as a parameter and iteration is carried out till the temperatures 
does not vary frem one iteration to the other. This final profile 
gives the values of and The velocity is known from 

inflow. 


2.2.2 Transient: 

I ■ mi I >i w ii i m wrn fl 

The side loss is not considered in the transient model, 
taking the lake to be at uniform temperature horizontally. The 
energy equation for the natural lake, without power plant, and for 
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the cooHDg pond, with a power plant, is obtained as discussed belows 
®) I^atural lake: 


The goveming one dimensional energy equatim is 


3 ^ _ f 

■^t “ Hih 2 

C7z 

with the conditions : 


C — 0} 

s 

at z 

= 0 

D 

at z 

= H 

■^T 

Q 


H^z 




= 0 , at the surface 


b) Cooling Pond : 


( 2 . 10 ) 


(2.10a) 

(2.10b) 

(2.10c) 


The goweming one dimensional energy equation is : 


^ ^ 3_T 


( 2 . 11 ) 


with the condixians: 


T = T 

s 

and 

II 

at 

1! 

O 

It 

and 

li 

at 

z = H 


(2.11a) 

(2.11b) 


w j (t. +a t) - tT + € 

} D s I 





at the surface 
.... (2,11c) 


2,2,3 ifetho d of So lution ; 

The cooling water body is divided into n horizontal ele- 
ments ( Pig. l), each of depth Az, The governing differential 
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equations are written in their finite difference form, see Appendix I. 
Separate finite difference equations are developed for the surface, 
bottom and intenaediate elements. These finite difference equations 
are solved by iterative technique. Iteraticn is continued until the 
convergence is obtained, indicated by the difference between T (l) 
and becomes smaller than a small specified value for all values 

of I considered, vihere r is the number of iteration and I the 

0 , s 

number ofy^element fron liie top. The initial guess T vl) was found to 
affect the convergence significantly, the closer the initial guess 
is to the final soluticn, the faster is the convergence, as expected. 
The following initial guess was employed, 

T°(l) = T , for all I (2*12) 

© 

For the transient solution, the analysis begins on January 1 st. , 
when the lake is in a fully mixed, or isothenaal, state and the solu- 
tion is carried out over the year. The iteration procedure is carried 
out till the temperatures remain unchanged from one year to the other. 

^•5 ^>^0 Dimensicnal Be circulation Mbdels ; 

A two dimensional model (Fig. I) for the study of recircu- 
lation assumes complete mixing in one direction. A vertical model 
wilh good lateral mixing is consideresd. This situation frequently 
arises in narrow lakes and rivers. The intake and outfall are con- 
sidered to be in the same vertical plane. The basic governing 




equations for tlie flow field, for constant fluid properties, eddy 
visoosily and diffusivity, and neglecting buoyancy effects, are 
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^ u S w ^ 

+ r— _ 0 

itjx az 


^ u tfU 

^az 


-iif 

^ ax mx ^^2 ^^2 


(2.13) 


(2.14) 


aw. a w 

u + w = 

ax az 


1 dp /B w a wx 

^ £,z mz J^^2 ^^2 


(2.15) 


ai 3 t /r \ 

•ii" ' + W i ^"' =r V. _ ^ “ "o + o' * 


u -r- + w ^y“ 

ax a® 


T +rT 

C C^z 


(2.16) 


where u is the velocity component in x-direction, w the velocity 
component in z-direction, 1^® viscosity in x-direction, 

}J^ the eddy viscosity in z-directicm, and Cg, the eddy diffu- 


sivity. 


The flow field can be solved independent of -Ole energy equar* 
tion. The equations are qiiite complex and require simplificatians 
before a solution may be attempted. 


2,3.1 ffifethod of Soluti on; 

It is assumed that the eddy viscosities in both the direc- 
tions are equal and given by, 

^ msi~ ^mz “ 

we may now eliminate the pressure termsfrom iquations (2.14) and (2.15), 
and the resulting single equation is obtained as, 



21 


-'..2 ^2 -^2 -^2 1^4 5 

“ 'i? ^ ' 'i? -i-a- ’ -^'az 




x" ^X 3 z^ 


) (2.17) 


A stream function is nov/ introduced j in. the above e<},uatioiij so 


that it satisfies the continuity equation, and is given by* 


(2.18) 


w = 

{?x 


(2.19) 


■Phna the Equation (2.1?) in terms of stream function y is given by, 


ay a (v^y) ^(vSo 


Bz ox 


ax3z 


= ij i^y) 


( 2 * 20 ) 


wliere . 


r.2i, 

^ “ i? * a" 


(2,20a) 


!Dhe above eqviation may be expressed in terms of dimensionless para- 


meters, particularly the Reynolds npmber, defined ass 


He s 


(2.20b) 


where, U is the average intai® velocity, d is the diameter or hei^t 
of the intalce Stream, and 1/ is eddy viscosity. 


OChus, the equation (2,2o) is expressed as, 

^ M ^ (' 7 ^) ( 2 . 21 ) 

3z 3x c>x -^z Re 


where, y is non-dimensicnalized by TJd, to obtain KjJ and x and z by H* 
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ITow Idle Equation (2*20) is solved for two extreme casea : 
a) When Reynolds number is very high ? the viscous terms may be 
neglected* The fluid then behaves as inviscid* For an inviscid 
flow the stream function satisfies the Inplace equation, vhich 
Is given by, 

s 0 (2.22) 


v/ith the following conditions, 

B'P 

x = o, o<a-<d , = n 

X = 0, H - d < z < H , ^ - u 

> S/L/ 

x=0, d<zi,H-.d, *0 

aw 

X = li, 0 < z < H , ~ = 0 
z = 0, 04 x 4 L , ^ = 0 

a 

Z = H, O^x^Ii,-^^ -0 


(2.22a) 


This gives the standard inviscid flow boundary conditions in terms of 
constant s’feream function on the boundaii-es, as given in Fig. I 3 . 


b) When Reynolds number is very small, i«e., the fluid behaves as 
highly viscoiis. This is true for very low velocities generally 
encountered in such flow, the problem being Imown as"creeping 
flow". Thus the viscous terms dominate and the inertia terms 
are neglected. ■ This gives Sqn. (2.20) as, 

vV « 


0 


(2.23) 
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wi1iL the following caaditions , 

2 = 0, 0^x:^L , s= 0, r*^=! a* Where, a is zero v^en 

Ox the surface i^ear force 


z=H, O^x^L , — ='^«0 
x=0, d4z4H-d ,|^*=|“ =0 
x=0, 0<z<d, “ = U 


Is absent* 


( 2 . 23 a) 


x=0, H-d<z<H, ~ « - U 

x=L, ,“=~=0 

* ' ^ ’ Sz 8 X 

The equations (2,22) and (2,23) are solved by finite diffe- 
rence methods, A rectangular mesh as shown in Pig, 1 is used for 
this purpose. At every point in the interior of the considered region, 
a finite difference relation is setup in accordance with the above 
equations. The various partial derivatives Involved are approximated 
by central difference approximations. Separate equations are written 
for the elements at the intalce and oCitfall.The surface el^ients will 
also require a separate equation, particularly for the energy equation. 
Thus the finite difference equatiois are solved by iterative techniques, 
-Thus, starting ftom the initial guess V ^ > we determine the succe- 

ssive iterates ^ . . V* • ...... H’. from the general f or- 

i,3 y * ^>3 

mulae given in Appendix I, Iteration is continued until the conver- 

— X'+l X 

gence is obtained, i.e,, difference between UJ . . and Mi . . becomes 

‘ .1,3 f 1,3 

smaller than the specified value for all values of i,3 considered. 


The following initial guess is used, 

s 0 for all (i,3) 

X 5 J 


(2.24) 



CHAHEE5R 3 


EESUIffiS AUD DISCUSS lOIT 

This chapter deals with the numerical results obtained from 
the thermal analysis of a cooling pend. The stratification of the 
cooling pond is largely dependent on the net surface heat exchanges 
■vdiich in turn is a function of the surface temperature. The net 
energy lost by the surface, Q, as a function of the surface tempera- 
ture and the wind •velocity, is computed by a numerical evaluation of 
■the various components of energy exdaange and is shown in lig* 2. 

The curve shows "that at low surface temperatures, the change in -wind 
velocity does not have much effect on the surface heat exchange, as 
compared to ■that at hi^ surface temperatures. At ■these hi^ values, 
a small change in the wind velocity, gives a considerable amount of 
increase in the heat lost at ■the surface. Energy loss increases with 
surface temperatures, as expected. 

The average surface equilibrium.: ■temperature, T^ , vhich is 
defined as ■the surface tempera'ture at 'rfa.loh the net heat exchange 
from the surface is zero, is computed fesr different sets of meteoro- 
logical parameters. In Pig. 2, "the interceprb of ■the curves on 'the 
x-axis, gives ■the value of -the equilibrium temperature for various 
conditions. The results, for varying ambient condi'ticns, are shorn 

■X 



in Pig. 3. Par eacsh variable, the remaining variables are kept 
constant at the values indicated as typical values in the figure. 

It is found that the surface equilibrium temperature, T , increases 
with an increase in the ambient temperature, as expected, and also with 
the relative humidity , ; since it reduces the evaporation, there- 
fore requiring a higher surface temperature for energy balance. It 

is also seen that T decreases with an increase in the wind velo- 

e 

city and in the cloudiness factor. This indicates that the net heat 
gain.'by the surface inc3?easeswith an increase in the ambient tempera- 
ture and a decrease in cloudiness factor, since clouds cutdown on the 
solar flux received. Also, the net heat loss from the surface increases 
with an increase in wind velocity, due to greater evaporation and con- 
vective loss, and a decrease in the relative humidity. It is also 
observed that the equilibrium temperature is more sensitive to ambient 
temperature than to the other parameters, such as wind velocity, 
relative humidity and the cloudiness factor etc. 

Pig. 4 shews the steady-state vertical temperature distri- 
butions in the cooling pond at various times during the year. It is 
observed that the vertical temperature gradient is slightly more in 
winter than in summer. This is jnainly due to an inc3:ease in the net 
surface heat gain during summer, Yhich is conducted and convected to 
lower layei:^ of the cooling water body. The steady-state temperatures 
are hi^er during summer months as expected. In this case, a unifom 
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tempera-uure in the surrounding fluid is assumed. The profiles for 

different values of h 5 the side loss heat transfer co-efficient, 

w 

are computed and plotted in Fig. 5 for summer. It is seen that as 

the value of h increases, the temperature gradient increases at 
w 

the top. Yfith increasing h^ ^ a greater percentage of the rejected 

heat is lost to the neighbouring fluid and a greater temperature 

variation with depth results. Similarly Fig* shows the steady- 

state temperature profiles for different values of ^ eddy 

diffusivity. It is also found that the profile becomes steeper with 

an increase in the value of , since this also results in greater 

n 

loss to the neighbouring fluid. These results were obtained with a 

constant value of h and ^ „ t over the entire depth. Hov/everj 

W il _ 

the temperature profile was also studied for h^ and C g ss func- 
tions of the depth, decreasing with an increase in depth* But these 
temperature profiles were found to be quite close to the previous 
temperature profiles for the variation studied. This steady-state 
model simply gives the temperature level at various times during the 
month and the general nature of ’the profile. A stratified profile 
in the surrounding fluid m^ also be used for more realistic results* 

The transient case is now considered. The time step-size 
is taken as one day and the ambient conditions were taken at their 
average values during this time Interval. Since the ambient condi- 
tions vary over the year^ the water body responds to heat rejection 
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from the power plant differently at different times during the year* 
The main problem is that the water body does not reacSi a steady ~state 
and its temperature profile continues to vary during the year. The 
transient behaviour of a natural lake , without heat re jecticn, for 
the measured Meteorological conditions at a location near Kanpur is 
shown in Figs* 7-9* Fig* 7 shows the transient vertical temperature 
distributions in the natural lake at various times during the year. 
Fig. 8 shows the variation of the bottom and surface temperatures 
during the year* The temperature profiles indicate that during early 
winter, as around the 1st* of January, 1he water body is in a fuHy 
mixed condition* V^ith the passage of time, the water body loses 
energy and its temperature decreases until the winter minimum tempe- 
ratiare is obtained* This occurs sometime in early February for these 
conditions* After this, the water body starts gaini n g energy and 
stratification occurs. Then, the surface temperature increases rapi- 
dly while the bottom temperature increases much more gradually, since 
the pond gains energy at the surface. This energy is gradually 
transferred to the bottom layers* The vertical temperature profiles 
in ISie cooliig pond become more and more steep at 1iie surface* This 
continues until the surface temperature attains the summer maximum 
temperature* After this, the cooling pond starts loosing heat* This 
gives rise to an isothermal top layer as seen in Fig* 7* surface 
temperature decreases rapidly due to this loss of energy, #iile the 
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bottom temperature increases slowly due to the gain in energy from 
the hotter upper layers. Dest 3 ratificaticia occurs in early winter 
and the water body becomes fully mixed, i.e., the isothermal state 
is reached. Ihe water bO(3y continues to lose heat beyond this point 
and its temperature decreases to a value close to that at the beginn- 
ing of the year and the cycle continues. The nature of profiles at 
various times is apparent in lig. 7 and the transient behaviour in 
Sig. 8* 

The occurrence of stratification is delayed by a few days, 
as shown in Pig. 8 and its effect on the bottom and surface taapera- 
tures, is studied. This relates to the delay due to the wind induced 
mixing of the water body. It is found that the surface temperature 
remains nearly unchanged, vdiile the bottom temperature increases 
with an increase in the delay period. A delay in the onset of stra- 
tification implies an equal distribution, of the net heat gain at the 
surface, vertically throughout the depth. Thus during that period, 
the surface temperature will be less than the surface temperature 
obtained without a delay and the bottom temperature will be l^ger* 
This affects the yearly temperature cycle of the lake and the bottom 
toaperature increases due to the greater input of energy into the 
bottom layers. 

How, the net heat exchange from the surface is computed 
by considering a sinusoidal variation of the surface equilibrium 
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■bempeiBture, , and using Eqn. (2.4), yearly surface and ^>ottaa 

temperature variation, for a natural lake, and the variation of the 
surface equilibrium temperature is plotted in Eig. 9, for two values 
of K, the surface heat exchaiage co-efficient. It is found that with 
an increase in the surface heat exchange co-efficient, the bottcm 
temperature decreases throu^out the year. This change in the hottcm 
temperature is maximum during early stratification and minimum nea 
the destratificatiott process. The surface temperature decreases dic- 
ing winter and increases during summer. Since the equilibrium tem- 
perature is less than the anrerage surface temperature during winter 
and is more than ihe average surface temperature during summer, the 
water body will lose heat in winter and gain heat m summer. This 
surface heat exchange win increase with an increase in surface heat 
exchange co-efficient. This causes the lake to gain more energy m 
summer, resulting in higher temperatures and to lose more energy in 
winter, resulting in lower temperatures. 

tte basic affect cf taat rejecticx wmM be tc tocreaae the 
aerfaoe tmperatuxe, *lbh «riea cror the year for the natural late. 
It will alao affect the hottcm temperature due to anaigy tranafer to 
the water body and toa totate water temparatuia to expactad to rtoe. 

10 ahows the aurfaoa and bottom tamparatura variation for the 
natural late aa wall aa for the oooltog pond, i.a., the late wito a 
power plant. Ihe diffarenoe between toeae two to matoly due to toe 
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rejecticsn, of waste heat fhom the power plant . to the cooling pond. 

The tempeirature level increases as expected. The effect of wind 
velocity, on the bottom and surface temperature profiles of the natu- 
ral lake is also studied and plotted in Fig. 10. It is found that 
with increase in virind velocity the temperature of the cooling water 
body decreases throughout the year at all points. It indicates thau 
the net heat loss from the surface increases wilh increase in wind 
velocity, as expected. The effect of the po-er plant capacity and 
of the eddy diffusivity on iiae cooling water intake temperature is 
also studied and plotted in Fi^. 11-1 2. • The effect of poi/rer plant 
capacity, i.e*, of the heat rejected, caa the intake temperature is 
shown in Fig. 11. This indicates that with an increase in heat rejec- 
tim, the intake temperature increases. This is mainly due to an 
increase in the heat, which imist be rejected by Ihe lake to the atmos- 
phere. It also gives the flow rate impact on the intake temperature , if 
the temperature rise in the condenser At is kept constant. The 
heat addition due to power plant is dependent upon w , the downward 
velocity in the cooling pmd and At, the temperature rise in the 
ccndenser, and a change in either gives a similar effect. Fig. 12 
shows 14ie effect of eddy diffusivity i or of increased turbulent mix- 
ing in the lake, csa the intake temperature. It is found that wxth 
an increase in the value of the eddy diffusivity, the Intake tempe- 
rature increases during winter and decreases during summer. Wiidi an 
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increase in eddy diffusivity the mixiiag in the water body increases 
and this increases the rate of heat transfer between the surface and 
the bottom layers of the water body. Thus a larger amount of the 
energy gained in the summer is transfejnred to the bottom layers 
which, therefore, reflects a higher temperature in winter, due to the 
transient lag. Mere energy is similarly lost by the bottom layers 
in winter and this is reflected as a lower ten^jerature in summer, 
for a higher value of eddy diffusivi'ty. 

The two dimensional vertical flow pattern in the cooling 
pond is also studied and plotted in Pigs. 13-14. Pig. 13 shows the 
stream line distribution for inviscid flow, vhile Pig. 14 shows the 
stream line distribution for creeping flow. It is found that in 
inviscid flow a greater portion of the cooling TOter body is affected 
by the flow as compared to that in creeping How where oaly a small 
portion of the water body is found to have flow at all. The part 
of the water body away frean the outfall of the hot condenser water 
is found to be undisturbed and thus in creeping flow analysis, Ihis 
part of the cooling pood does not play a significant role in heat 
rejection. The flow pattern, in both the cases is more or less as 
expected. Per a very big lake, the area where the flow exists must 
be determined for a proper understanding of the flow mechanisms'. 

The problem may now be solved for the temperature field. It is ex- 
pected that liie temperature will decrease vertically downwards and 
away from ttie outfall. The adiabatic condition is to be employed at 
the walls and the surface energy exchange at the surface* 


CENTRAL L IBRAWf 




CHAPTER 4 


COHC333SIOIIS 


Pram the discussions in Chapter 3> it is noted that the 
results obtained are physically reasonable* The steady state pro- 
files show that the water body is nearly stratified throu^out the 
year, v/hile the transient profiles indicate that the cooling water 
body remains in fully mixed conditiaa during early winter and in 
stratified condition during rest of the year. This is mainly due 
to the time dependent nature of the water body* The transient pro- 
files obtained by taking a sinusoidal variation of the surface equi- 
librium temperature are found to be similar to those obtained by 
considering the surface equilibrium temperature as a function of the 
exact net heat exchange from the surface. The equilibrium tempera- 
ture in the latter case is defined as the surface temperature at 
which the net heat exchange from the surface, at the given ambient 
conditions, is zero* The effect of the heat rejection from a power 
plant is found to be a rise in the intake temperature, #iich shows a 
greater change for a higher power plant capacity, as expected. The 
effect is found to be larger in summer, due to the higher ambient 
temperature, vdaich nnkes the additional heat less more difficult* 
The effect of a variatiaa in power plant capacily, in turbulent 
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mixing j as reflected in the eddy dif fusivity y and in ambient condi- 
tions is studied* The results obtained prov^ide guidelines for limit- 
ing power plant capacity for a main water body and all indicated 

the expected recirculation effect on intake temperature* 

* 

The two-dimensional flow pattern 9 in the cooling pond has 
also been studied for the two limiting cases of inviscid and creep- 
ing flow* The foxmer applies for large outfall velocities and the 
latter for small values* The present work can also be extended for 
other flows that may arise in practice* The temperature distribu- 
tion^ #iich is closely associated with the flow pattern may also be 
considered for future work* 
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APEBUDH I 


OTIIE DIPPBEENCE EQUAIIOIJS 

The finite difference equaticns for one and fwo dimensional 
models are developed as given below: 

One Dimensional Model ; 

The equations are developed for surface ,bott cm and inter- 
mediate elements separately, with the follOTdng parameters, 



ft 

(1.1) 

B 

' , „ 

^ rfb. a 

(1.2) 

0 

h 

w 

* PC L 

i P 

(1.3) 

D 

1 

“QO 
^ P 

(r.4) 

Hovr the Eqn. 

(2.9) is represented, in finite 

difference form, as : 


a) Surface elonent. 

w (T (n) + at) + B . T ( 2) + C.T - QJ3 


T = T (1) 
s 


w + B + C 


(l.'S) 


b) Intermediate elements. 

(w + S) . T (I - 1) + 3.T (I + 1) + O.T 

6 


T (I) 


w + 2.B + 0 


(1.6) 
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c) Bottan element. 


(w + b).!D (n - 1) + C. T 
I. = T (n) = 


b w + B + 0 

Eqn. (2,10) is represented as : 


a) Surface element. 


T = T (1) = 
s ' 


B . I (2) + 0 . T (1) - B , Q 

B~rc 


(1.7) 


(1.8) 


b) Intermediate elements, 

B . T ( I - 1) + B . I (I + 1) + C . T (l) 

Id), d.9) 

c) Bottcm element. 


I - T (n^ - B , I ( n - 1) + 0 , I (n) 
b“ B+C 


(i,io) 


and Bqn, (2,1l) is represented ass 
a) Surface element. 


( I Xn) -t- at) + B , I (2W 0 , g (1) « 
b) Intermediate elements. 


I = T (1) » — , ^ „ 

s ' w + B + 0 


m ( t \ , Iw B) , I (1 ~ 1) + B,I (I + 1) ^- 0,111) 
^ " w + 2,B + C 


c) Bottcsm element. 


* T (n) = 


(1.12) 


(w -t- b) . 03 (a - 1) + 0 ,0! (n) 


w + B + 0 




(1.13) 
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Two Dimensional Becirculation Llodel : 

The following non-diraensional parameters are used, 

* u . d 


(r.i4) 


— X 


z 


(1.15) 


V 


i,3 


2 + 


Now the B^n. (2,22) is represented as , 

. , V ,. ■!.., 1 ?3+ T i i3z.lL 

• O 

^ ... (i.'te) 

( 1 . 17 ) 

and Eqn* (2,23) is represented as, 


where == (V^) 


VU (HB.AA - DO . AB > 2. Cl . IB . AO -f BC » AD 0C.MJffi) 

^i.3 


(I.I 8 ) 


*ere, A* - ,j + Vl-1 , j 


“ = Vi+2,3 * Vi.2,3 


AC 


“ V i+1 ,3+1 ^ V i+1 , j-1 V i-1 , j+1 V i-1 ,3~1 


AD = 4 ./ . . - + W. . , 

' 1,0+1 ' i,0~1 

AS ss W . . « + 4 } . . « 
• 1,0+2 r 1,3-2 

« (L/H)^ 

cc .= (L/a)^ 

X « ( D/E) • A X 
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m. rr 

( A 



IB = 

( . 

Az)^ 


DC = 

(Ai)^ 



BB = 

4 . d6 + 

4 . 0^ 

. IB 

li 

o 

4 . DB • 

C^ + 4 

. CC . Di 

BD = 

6 . DC + 

o 

• 

CO 

.336 + 6 


(1.19) 




<1* t ; t" i 


f; ^ h^'^i *; %> ^ i'lr^ i p, u k ^ u U fe >.v 

#< '’i ^ "'ii'f ’ ,^. 7 '*, ' ^f»,’''i'i" , '; ;’i ' v'i' 

!-> r 5 p *#': 'I,cnrcr^ ', M'“ 'pii ',m - ”^>‘.%t'U'^ i j' ^ f/, y vt,y !V‘- I a ^ : 


a f%' 4i., 1^ ^ ^ ^ .« ,, .^'y y 'a^ .i^s fe- f/ ,ig * ^ ^ ^ 


'y # 'f ^ 4 |t. %* -is' ?<•; f '||;^ ^ M 'r# If 


-p it ^%y< ! yi/ $ypr‘ ^yyyyy o^*" 

































































MWMWjllB! 

n^HBi 


f'r/fi^irf ,1^/#^?^' #44t^ ' M" *«^^i’^^.i uv/t«' f-=n^ 






























>» ^-^*4# *-»^‘ t'^*'* * ’'# ^ ?i" -i, |‘^f> f- v^ » r * H 


>? ^ f^ '*^ (*^»i, "" "“ *^f"i''*’^^y &■ ■""rif* # ^ ^ *>' If I ^ f» ^4 4 * ^f> ^ ^ n U ^ ^ A^'1 

» "m .hi^, 1^! i ##;>ll4r4V#^W#4/.Tis4i^m‘i^jlr L '<4^? " JteS- 1‘^ 









* ij. . ^ ^ 


::i Vs fs ^ t ; -v^'^ ' rf h£ f; . y v%^c^ . « v uA7^ui * 
















'■llM'<^#'V ^ i*.,-i‘/-:'t'^#''*'''^^^ a 5 jIs; r-: '''ff - ;fs;' *' ', ,; , , 





























' t Y t 

"V ’* *i i'^/ 

# ^ 



^ ^ ft 04 

I^^P 
























iii®»l 





























j1 Y 1 1 slfl/ 1 ^1v '|‘ "'^,1'';' f y"F 1'" '''^1't1v''^''^''1 'f f 'If if 11""f^'';f' /""'j 




























,-,v., ..^r 




V# 













f .. / % f ^ .w ') 4 '^ ^ * 









*1^'' P''" '* >' ^ ’ f|>/f ' I . 




h‘',;P , r^'^f 








>V ^ ^i# t *'4* (t r^ » fe Ai & , fi ^ -i i fe . ^ ^ X* ^ frij ^ ^ ^ 














Ij') «f=a^ >^5-‘<is * '« 'ffcl'^ m ■{ * ft J- & *' ^ rf V« I?iv fSi ^ i!*, S(|- ^ ® >C t< ® 'f'V iU ^iis ^ iiiS s&'^^S ^ 

ff'' ^ 















^,., 









»ff,' •:.. ,- •'. •• '."’’‘"S' -ft V', w-a . - «• -fc (. ♦'i j. ->'‘i„***> -j, * *.•*..». , .' j.' .'. ST ^ 'f" j) 

























If,. ^ 










S' *'cv'';^-i>vic//, <;■■>> V': ',? r :a t;?' 

V4< i'''' • %■ ''»■ ,1'#^--C'^'i»’ vt^ 

— \ -fir' p' ^ ^ ^ ^ ^ ^ if t 5 5^4. «. « Si , 















• ,•» f#.#S;/'-’‘4!^ lilir t#'" 4V*i>y;?i€||J-f,, 

M 

mm 

r- , f ^ s 

14^ #®'- m ^ ^ V f "i V « *f g^ f ^ v},# f >*. g ^,jt 

f > ^ ^v.v’ri^ ,^44 . U r4;\ ? 






" *» ^ »# m* ^ '-T ^ ^ •vT '.ir «r 4.f gi iff* ^ -^r ^ <*$" “v?- ^ vv ^ 4' 


;* « 






' - i«;' ’-f '■' '*■ 

■ ^ ‘ -fit' ^ 

«« *. ";.. i ••4, 
















^ ^ * .. *. .. .. 4^ *. ^ M ,4 .. ^ .y ^ 4r,. 


-.7’ 

W.f^rU. 
















rX '^1^77 , ,,--7 ,>1 y^'f>/,wV^>^ . ' 4t^ tt ; .r7; 1-’^ 

' ^ , ‘7"-* ■"'" ^ h'^ /^-r, 4sv A^if , ‘T'l'- 'fv 
















;. tf 7 ^ \, -1^ )^ ff t/ f 4, '^;tV,^ ^ •‘7 














^ , 

















































l^■!■ 




nif^i, , » >i y V 



// ip .'-V4's,| ?«■ 

>>3|, *% '# § , s ^ ^'i'iiim^ ^)> f S' p4^#-^ S'" j ■'*«'“ % yfi'^'ij j# 

■ ?;■ •:■''■>?■ .'«4 ri! % ‘t^ri# :y;m 



‘ : rh,. ^ ^ i/l4 # i> «f 

! '' '# s ’'/’I" ' . " ? 1- / t \i4't ' ‘ ^‘"1 / " . '-''^ 4 ^ y' - %fi‘l 




:■ £i f : :■ « 


*• '1 vi>' ■ ' .',Lr-, 




jif^ i^^'-f-vi ri V ■w V .<f ■8i),;i’ •' 0 ' >fiypir- ^ ^ 

, 'fi ?j is ^ ’% |y ^ V «.,• <'i|' / 

/Ji^ r r^-^: ^ ‘ ^ PP,"'? " ' ’ ‘if. / .W% l,Zni\ 

/4 4 , f i\,:\ \ ‘^y ,<,f >ft 






" ’# f iflJ-tP: 


^ ^ ^ 4 T ' 1 # ' ^"7^. 

i’^Jfw."^' *- ^;/ ' .V, ' ,^,.-( * V'#5j;''''VPa';ft,-;!*’’ Vit,r<j r ■.„;w-?*'’i- ',# 1;;pJ 

•' ' -f.p*'''^ "".‘ #^'T"'"''' ■ iW'"’ •■tt;.-'!.'/.' '•k"'>' f"/! •'• .',t', 'i,"^'/,J.'y- #■ ‘,p.5 ': 

t^^tf'l ^ i' *.v ‘ l^'-.^Vi? ; >’ ' .^'^”'-'V “, 1 ' ', 1 ' 

■ __.«. (*, - ■tfi->«J.--^)f,,;,i«,«,«(Jl(.f>»t. 4 -«»<SlffeD> 

^ 










“ A; ’7 '>f''^./^ /aP/7:p 7 7 7/f ,i^; pM/fpp.pp*-' ,C7ff . r'c^v, 


Pr|/.t7 'r^j 









































•? fe 


< s ^'- 


‘^’ <^rA^ ^ '" '^ ’^ 4" ^ «-*^ ^1 ? Hi* # ^ V '#.#> /, ^ ^-|:r 



' 


if .« ii/ u ‘.tPvf'^ r'’,\l^.i!^^?i^',:i^ !iku^'*^”^'^ '‘!4y^' ft 4H V ^/<Hf , j ; 








'i^'L ^ 












li^SPiWiw 

, 

1^™* 


^ u ^1 f .i, * «. >,4' '^ ”;4 a fv '% iv « 5 |. f' ‘^‘/u ^' ';v5f ,' 

' <?* ^-f:, /> ^W'' 1 */ Y'^> ?^ - V *^ii‘ V >ti'’V<t Vsjt''^|«'^P;/,v hV *^if^'4iH\‘f^|- 




















^Ps-?./ P" ^ '^ ^" .ivjvp%^y. /'?'^. .1 h'^ / >itiV ; /\4/^' //V 



I''?::?"'*'"-' >^,'1 F. 




















■*h-t^i V -L-^j < I'" 









# 









; l)>f -i-y. ♦:>^rfe^^:^: \v ■ ^ 1 >^ , ^ . 












v.-f ,i* ^ ^'v ifr*',/ 1 ^ 

<*-^ x'‘f V sr «r # j , ^ ^ ! 




^ V ^p -ip 



4.- i„ 5 ^ #■' %T '4" ''^ i -ff 




^^^^ xmWim XiM^: 


vr,3-''’';,t >?'"■/ T ’-y--'r ^r-T^r •;?* -cs'-ff 'i"--t <'« -r •’* -:» *; 













ift»8f^ttipj4|#iS8ijS 

?’'5? » ,i' ^ r ,1 

inv 



[,f'^' I'Vi 





4/(f if\Zl^ . 









i» j f"^ 



^ ) "” A / ' 







*/* ■f'.w ' 




pi®w 




r4 ^ ^ ’# 






• . pyxd . 








“iV --^ 4 i^i ^ 4 ‘^ t" ^ V ^ ^ ig VT^ w ^ /i^'' ^ m «. yi ' 11“ 

4^ f«-vt f^ f^^'vt. ^(^|1«s,.^‘s^ |*»^^ '-w,|tvt| «*,|Si‘^'" , 






ij>. ’•■>'4, f% ^'S^|!*,*|%-«^*^.t|<>^»i "'^^ 

5^“ tiij’?!-.' vi.“> '*^"nv :,r* ‘»a X v* if^x^x^'s^ # . 

>$1 4' ^ 4 # fe^ ^5^ ^f. ■*" 1c. 'i^ t, (fc* f ^ 4, f? ^ lirf /f p ‘ 


IT ' , "7: r'?/ ;x-" "" «"" ^ 

^ ‘ mI> ; -,r ^ ^ # 5. ^ i** <,» ^ ^ s/ ii/ ^ !gr itV 

Tr X-^' ^ *;f»-’‘“‘# .-f -t - ' - t ?'«/"''4? *>TVu.f fsr' T^ “"'-f </ > 


;.* y ti .^X’ ^ *’ip'' '.> 4 r -X. "X'. ' 4 < -r^- f y^X' '4. /t"'" 4 

























iOooooortor>r^r»niot^r^rjnior>ooonor«r>rif*»nnono 


APPS-4DIX III 




OETe^HINATIOK OF ONfc 0I«Eh!S10r4AL VK^liCAL Te!''»FERAfUr.c Dl STPlIBUTI 0 
IM CCOtING PCNO 


!»««*»#»» iw® «>« «i 


THE INTAKE AND CUTFALL ARE COHSiOEREO IH SAMP VERflCAL PtAHP 
TH| INTAKE IS TAKEN AT THE BOTTOM OF THE COOLING PO'^-O HHILE THE 
OyfFAiLIS TAKEN AT THE TOP SLRFACI: 

A NATURAL LAKE IS A COOLING POND WITHOUT POWER PLANT 

TOSUL IS SUBROUTINE TO DETKRHINi* TOTAL SURFACE HEAT LOSS, FROM 

COOLING POND 

¥IfiHA is SUBROUTINE TCI OETERMiNE SURFACE EQUlilBRlUH: TeHFERATUR^ 

OFiZ IS DEPTH OF ELEMENT CON SIOSREiD IN VERTICAi OOWEiW.ARO DIRECTI 

DELTIM IS TIME STEP SIZE TAKEN AS OME DAY. .. 

re IS AVERAGE SURFACE EQUILIBRIUM TEMPERATURE 

AHUH IS F0DV THERMAL DIFFUSIVITY 

SETA IS ATMOSPHERIC RADIATION FACTOR 

CF IS CLOUDINESS FACTOR 

Y IS WIND VELOCITY ALONG WATER SURFACE 

QTS IS TOTAL SURFACi' HEAT LOSS 

OSF IS NET SOLAR RADIATION ON EARTH SUEFACG AT -KANPUR 




tn <![Mt »«1» «IM iMt ^ 


STEADY STATE -SOL UT I OH 




C MONTHLY TEMPERATURE OISTRI8UTIOM IS COMPUTED 
C . . XOOLIHG POND WITH POWER PLANT IS CO'^iSIDEREO 

fIBJQS 
$IBFTC MAIN 

DIMENSION ATA(12ltARH(l2I,AVfi2l,ACFa2I.AB£TA{i2J tAG$F.C12> 
dIkInSION AESTAtlll 
OlMENSION Tfl25|*TOLDU25J 
COMHON AtA*APHtAV,ACFjA8ETAtAQSE*A|STA 

R|A0 TUIatAC i! },ARHCII l,AVt IU,ACF| IIJtABETAf 11 } tACSEO 1 1 , AESTA { I 

Ni»N*I 

RO-IGOO* 

OElZ«|*S 

CP*i« 

00 22?. 

TA^ATACltl 

CAiL VISHAlTEtliJ 

AT«T| 

TB»ti 

FOK-TE 

oox»Te 

C« (HW’e'BELll/ltO^CPm I 

X2*M4a 

X3»W42.*B4C 

^ It«wf 

do. ft ' 


fl»*i 


©f LT l+BtFOX+CW-B^OTS >/ X 




to 10 





Wf 


00 00 I ® . 

D | FF « ABS ', . 

IFtBIG . LE.i 
COHfIHUE 
IF(8I0.te*0,0l) 60 TO 4Ci 
IFClia.EQ.SOOJ GO TO 40 


cii-Totoc nj 
OIFFklG=OIFF 


53 

30 


ff 

. 6.0 

18 

• fl' 


GO TO 
40 PRINT sSjITR 
PRINT | 9tTE 

forhatTsN . s ) 

■ FigM4Tl2X,*TE«=*f F6.2J 
50 FORWATCr" 


PRlNt 

fqrmIt 

FP««T 
FORMAT 
■ STOP- 
END 
C 

$I8FTC TOSOL 

SU8R0UTINE 



SUBROUTINE TOSUtieTS.Tl, QBR,QS,.«H,Q", I U 
DIMENSION ATATl2}.ARH{l2J,Ay(l2l,ACFCl2l,A8P:TA{l2J»^'0SF{i,2} 
DIMENSION AESTAdlj 

tat- iktutnf mk .Mtt * MMt '« ii<w if* d» ^ itvi. « ' - 


10 


2 ?> 

30 


n 

92 


50 

40 


60 


DIMENSION TT{5IlfESTTC5lUP01t55J 
COMMON ATA»APH»AVf ACFtABETA, AOSe,AESTA 
TA*ATA( |I) 

ESTA«A|STA(IIJ 

SSE»AC3$£tIll 

8EtA«AB|TA(Ii> 

CF»AGFf|lJ 

V*AV(li1 

ffTtHL*’ 

DATA ESTT/%204, 9*839, W.Sl2,il»225 
1.47, |6#47, 17*53, 18*64, 19* 82,a.&6, ‘ 
a.53,||*8l,33.69, 15.65,37*72,“' 
3 * 31 , 50 * 33,6 
4. 54, .102* 10,.^ 

DO 29 I*l,5_ 

IFiTLGT.TTrm GO TO 29 
IWNT»! 

GO TO 30 
CONTINUE 

IFdWNt.fG.ll 60 TO 91 

42 » 1NNT+1 
60 TO 92 
ji*IWNT 
'2*lWNT+2 
0 40 ■4«41,42 


A *3 "k ^ . A A e ie: #4 ^ 


,15: 

. ^ ^ , ,3!3 

69,52*44, 55 



IcLfji - tPiuSfSdl - TTnHl / tTTfJI-TTdn 

CONTINUE 
EST1»0.I 





♦ESTT(] 
'73.+Ti: 
s #* 2 |; 


i 


RETURN 
END 


IIFTC 



#*4-BETA* f 2? 3* -HA I #»41 

*llll-RH#€STAl 
"" 'TAl ' 


,21 , AV c 1 2 1 , ACr 1 12 > ,AieTA t ttl # i 


TSt 

tUlsPHl r, TF C 5 1 H ep 
V, ACF, 


)fM 

CF=A€ 



os|=Agsg<iii 

ESTA«&ESW(UJ 

TSCiJ-W. 

DO 2S I=2i5i 

2S TS ( I J*TS r J *^1# 

DATA 6STS/9*20A*9.839, ly,*5i2, 11. 225t 11.982, 12.781, 13.628,14^524,3.5 
1.47,16.47, 17. 53, 18. 64, 19.82, 21. -'se, 22.37, 23. 73,25.20, 26. 28.34, 3^ 

2.63, 31.81, 33. 69, 35. 65, 37, 72, 39# 89, 4?. 17, 44. 55, 47. 06, 49. 69, 52*44, 95 
3«-3|, 58. 33, 61 . 49, 64.02, 68.25,71. 87, 75. 64, 79. 59, 83. 72, 88. y'2, 92. 52,97 
4.54,l02.W,il7.2l,li2.53, 118,96,123.83,129.85,136.11,142,62,149.4/ 
no 33 1=1,51 

Ti*fscn 

FSTi»ESTS{IJ 

'QBR»4.74E»08tC C 273. +TI 1**4-6^14*1273. +TA}#*4) 
CS«-Ci.-C.eO?l*CCF**2}.J*C*SE 
Qfl«9.7*f0.48+6.II755*¥)*CT1»TA? 
Q|*l9.4*(0.48+0.075fi*VI#C£-STl-RH*ESTAJ 
OTS«QBR+CS+C-H+QE • 

TR<1 }=T1 , ' 

STSRCI)=QTS ■ ■ 

33 CONTINUE ■ ' . 

DO ICO 1*1,51 

IFCQTSROJ.LT.O.arGO TO ioa 

IWNT*I 

GO TO 101 

100 CONTINUE 

101 Jl*IWNT-2 
a2*INHT 

DO 102 I*J1, J2 

ggto:Si,j2 

CCNTIHUE ■ ■ 

TE«0.O 

DO 104 K«JI|42 ■ 

TE*TE-i-Pai«Kl*TR{KI 
RETURN 
END 

fENTRY 

§ •" 


103 

102 


104 


g 


TRANSIENT SOtUTIOM 


:iitS;atc#;j^S4*4i|*4!4*4**3|<9|t**#***********4*i|!***i! ******************** *******#****'' 

OMiY NATURAL LAKE IS CONSIOEREO 

ittg^ 1A5!r .*#|ps *^^-'^** *^p!» •W^P* WlP •»*«p 

118408 ' • ■ 

fllFTC MAIN 

DIHlNIjOH ATA(3651|ARHC365j,A¥(3651tACF{365HA8ETA(3651,aQSE(36S3 
lOf^ , AE.STA v365 ? , T E R.S T C I T 


1-5 1, -f fig] 

wunnON"AtA,ARft,AV,ACf=iA8EtA,‘AOle,AgSTA ^ ' 

REAO 100,fATACll,lRHl!l,AV*! ItACFCII ,A8ETACI i ,AQSE f U ,AERTAC 11 , I »1 

H¥H « 

mk% 

M ||p4.0 
Ri*lw0.s 
ITR*0 
CP»1#O0 



I ™ .'i * i'l 1 
= C8«ni 


^=T C^!lHC=i=T(^'i S/Xl 


7 

7! 


80 


%T~T ( ■! J 
”AVRC^--=T(iJ 
70 3 
TCO 
T3=( S'-. 

M J =TB 
JK=1 

nn 4 ir=?, Si 

IF(T{il“*r( IH6, 6,4 
J K=I 

TAVRG^TAVRG+TUK) 

C3NTirrJE 

GO TO 9'i 
T{IJ=T4VRG/FiOATCJKl 
DO 7 1-2, JK 
T{n=Tii! 

CONTINUE 
1 j 

i FCC-;;r8l?*LT,<J.,riI GQ TO 
PRINT 


-i)+BOTCi+ii+c-;=TCin/K 


IF{J,NE«lf Gl 


IT 


TO 


iO 


tCCiCUl GO TO 99 


99 


iCj2 

10© 

101 

77 

41 


=ITR+1 

-aG-0.0 

00, 1C 1 = 1, N 

CNECK=ABS(Tt i J~TFRSt'(i H 
I f {DIG. LT, CHECK J BlC=CHgCK 
CONTINUE 
IF{BIG*LT,C 
DO 5 1=1, N 
TFRSTd J=TCI I 
CONTINUc 

IFICOOE.FQ.l. J GC T^i 98 
IFCCODE.GT.Z.cn GO TO 77 
CQCg=l0.f5 

PRINT 101, J, tTCn, 1=^1, Nl 
0 C® 3 
60 TO 98 

FOPHATflH ,7F9.3J 
F0RpT{?F9.3l 
FORMAT (IH ,I3,15F7.3J 
PRpT 41,1TR 
F0RMAT{1X,*ITR«=I=,IS) 

STOP . 

END 

C , ■ 

SIBFTC TOSUl 

SUBROUTINE- TCSUL(Q|S,Tlj,. . 

OIHiNSION TTf51J,€STTC5i),P0ll55l 
OIMENSIQN AESTl{365^,TFRSTa5),T{ 151 

DINENSIqn ATA{365),ARH(365},AV{365),ACF{365J,ABeTA(3l35),AaSE(365) 
COMMON ATA,.ARH,AV,ACF,ABeTA,AQSe,AESTA 
P.H®ARH(IIJ 
TA*ATA(HI , 

V®AV(1U 
CF®ACF{II| 

9ETA»ABFTA{I1) 
os|»Aase{iu 

.^■SfAsAESTAtlll 
T-TCl)*!?). 

DO 10 1=2*51 
' TT( II»TTI i-U-frl*. 

DATA ESTT/9. 204, 9. 839, 13*512, 
l.47,16*47,n*S3, ■ 

3.3|,5|.33i6|-49,64,8i,68.25,7|,87,75.64|79.59,f3.72|88*02||2.f2, 97 
4« 54, 16 S * ly , 107 •2lfl,i2#53* 11& *56* 123* 83,429*85*136* 1-1 , 142*62 , 149* 4/ 
00 16 1*1,51 

- IFif|*GT*7T{IJ) 60 TO 20 

GO TO 30 
CONTINUE 

IP{IWNt.EQ*tl'-«Q TO 91 

3i* 


in 


IT 


ZQ 

30 


91 

92 

50 



30 IS- 4*31, J2 

31, ( 4 )* 1* 

^ 40 



TKIHI/dTiai-TTUH 



4r' cirvTifiin- 
ST '‘a-’- 

60 •;:; ST ST i * p t ’l ( i } f st t { i j 

() m =4. 74 H". •- ( C 273 t +T T 5 »4- T'liTAO { 'P . « + rA H «'4 1 

C S=- C i *“fl • 1!) . 7 1« C C F 2 } > S -'V 2 4^ 

. "H=q.75»'{ ■i«45;,40,'.(753 = 'VI=»‘{Tv-.*^AI • 

• 1 9 « 4* { 0 • 4 8+Q * 17 5 0=^ V ) 1= ST .',• S r A ) 

QTS=17.’^+eS + CH+C!r 
RTTURN 
PNH 

$Ri\'’RY 

C 

C COOLING POND tilTH POWER PLANT IS COMSIOfcRED 
SI8J0B 

$iBFTc nm-i ’ ^ 

DIMEIiSJO’-^ OF ITT C 101 
DirEf^SIiTl A?^UHHC ITT 

n I P.ENSIO.’! ATAC 3651 * ARH { 365 J , AV{ 365 ) , kCF t 365 > »A8FTA { 365) * ASS'^:: {365 I 
01H,”NSI0N AT.STAC 365ltrFRST{15l,T{ 15) 

CUf'MCN ATA*ARHtAVtACF, ABCTA, AGSE, ATSi'A 

Rc’AO ill, (AT4{ I ),ARH{| ),AV{I l,AGf Ci) ,aB-TA { 1 1 , AQ$f ( U t At .-.TA { 1) ,I “I 
i,365) 
tE=15»8 
AMUH=0,036 
DEL2=2,0 
CyDF®l*u 
DSLTIM®£4*0 
RO®lC0i1*O 
!TR»0 

cp«i,oe 

0ELT®9,« 

w=o.oii 

C®DeiZ/DELTIM . 

X3®I*/{R0=«‘CPS 

JC»£ 

■ N1*N-1 
3=AMUH/DFLZ 
Xl®W+B+C 
X2«W42«^B+C 
■ DO 1 1=1,15 
TFRST(I)«TE 
T{|)«TE 

1 COMTINUE , 

AT«T{1) 

TB=T{«) 

98 DO 2 J=4C,365 

T(i^=ISI|ySDiltf]4i*Tf2l+C*T{l)"-QTS^X3)/Kl 

fkl~ini 

TAVRG=T{1)- 

DC 3 1*2, N1 

3 . Tfl)*(«W+8l=S‘T(I“*l)4a*T{l4U+C’»=TCI))/X2 

TB={ {W4B)4T(N1)+C««TCMH/XI 

T(NJ«TB 

JK=i 

DO 4 1*2, N : 

IF(T{i)--Tf I) )6,6,4 

6 

TAVPG=TAVRG4T( JK) 

4 CONTIWe 
iFtJK.Eq.ll GO TO 97 
TCl)*TAVRG/FLOAT{JK) 

DC 7 I»2|JK 

r cn«T{ii 
7 CONTIf^UE 

f- TsiT ( 1 ) 

97 IF{CQ0|*L7*9,S) GO^TO SCI 

8u iFl4«rw£#3*i 60 TO 2 

ITR»iTR+l 
B16=5*t 
DO 16 l*|tN 

CHECK»A8| CT { I l-TFRST{ I ) ) 

10 

IF{8t6.LT*O.O0OlI SO TO 99 


5 



Cl a 


Jli. 'Sil 

b9 

i; i 

7? 




IF{CPOjl«KQ*.l. 
IF(CCa^*6T«2. 
CODE=J;;jr^‘ 


) GO T 
GO TO 


*7i?- 

' 77 


01, J, (Td 


I'M' 

J C~2 
GO TC 

FCriWATClH ,7F9,3) 

FOR’VAT(7Fg«3i 

rGPJlAT{3F7.A) 

FORMAT CIH ,13,lSF7.:j| 

PRTMT 41, ITR 

t-QR‘'?'A’^dX,=S‘ITR=«^,I5f 

CUNTI.'IUli 

STOP 

F;MD 


liSr-TC T 

c 


;TA{365I ,a,0SEC36§J 


IQ 


SUL 

SUf'RCUTlNF TLSUKQTSfTl.in 
OH^r-NSlON TT(51l,l”STfdi)rPnL( S9) 

DIHEIiSIOf! Af:STA{ 3655,TFKStCl5l,T(15J 
D I Me :v: S I ‘J-’ AT A ( 36 5 } , ;• RH C 36 5 J , 6V C i, to5 J , mCF C ;-j65 1 , 4BI 
COMMC^ ATA,ARH,AV,ACF, A8ETA,AGSE*AESTA 
TA=ATA{II) 

V=AV{II} 

!"STA=A6STA(II) 

QS£=AOSHiII} 

BETAsABFTAt in 

CF=ACF{II) • ■ 

RH=sARH(IIJ 

TTUJ»1Q* ■ 

DG 10 1=2,51 
T|Cn=TT(l-"ll + l. 

0 ATA, eSTT/g. 2Q4, 9. 8:.9, 10, 512, 11*225, Il»982,l2» 181, 13.628, 14, 524, :5 
1.47,l6,47,17.53,lft,64,l9*82,h.-.-6t22,37,23.75,25«2S,26,73,28.34,3j 
2* 03, 31*81,33,69, 35*65, 37*72, 39*B9, 42*17, 44*5f?, 47. 06, 49*69, 52* 44, 55 
3. 31, 58*33, 61*49, 64*f'0, 68- 25, 71. 87, 75*64, 79*59, 83. 72, 88, 02, 92, 52, 97 
' :i, 112*53, 118*06,123.83,: 29.95,136.11,142*62,149* 4/ 


GO TO 2Cf 


20 

30 


r. 54, 102.10, It;?., 
no 20 1=1,51 
IFtTl.GT.Tttin 
IWNT=I 
GO TO, 30 
CONTINUE 

IFdWNT.EQ.n GO TO 91 ■ 

4l=IWNT-'l 
J2=IWNT+1 
GO TO 92 
91 J|=I«NT • 

J2=IHNT+2 
92 00 40 J=J1,42 

00 40 I=J1,4 2 
IFII“415e,40,S0 

SO m(ai*cpoLi4)^m«^TT(n jj/fTT( Ji-Tiun 
40 continue 
:: EST1»0.0 ■ 

Ei?i=Esll+JSL?IJ*ESTTCn 

CiER«4.74e“€!8*((273*+Tn*44-BETA'«>!C273. + rAl’i‘»4J 
QS=~{ 1.-0. 0071^ {CF^42n*QSG/24« 
a£=19.4*(0.48*i-0.0755*VI4« ESTl-RHl^FSTA) 
0H«9.7*e9.48'»-0.0755*VJ*Cn-TAJ 
, QTS,«CIBR+QS+QH+QE. 

RETURN 

END 

SENTRY 
C 
C 
C 
C 


RERIOOIC SOLUTION 


NATURAL LAKE WITH PERIODIC VARIATION OF EQUILIBRIUM TEMPERATURE ' 
IS CCNSIOERED 

BJOB 

OSMINSION TEASIMSI , '' ■' 

DIKENSION A|StAC36SJ,TFRSTn5ltTtlf» . ■ ■.-■ ' ■ 

OlMENSlGH At A{365J*AftH(36S»,Av!365J,ACP«3«SJ.,A8itA«36Sl,AQSE(365J 



ATJ», ABETAf ACSI;;,ftFSfA 

3iH.-:r^’‘naN aissi •>65 5 

DiM^r'SIcrj T 

' e'AO l.y C ATAC I ) sAj^HC 1 1 ,AVC n f ACf UJ ,A8.. 1 a ( I } , AQSA{ I J f i 
If 365 S 

nHiz=r. 1' 
c-nc»i=i* ■• 

;H;LTIMa24.»^ 

^OslOOO* < 

ITP= ■ 

CP^ie^O 
AS«lUH=0,-''’36 
B = AHUH/0£;'lZ 
C-D£LZ/OELTI!^ 

xi*e-+c 

XZxZ.’S-B+C 
^3aI«/(RO>*CPJ 
JC = l 

no 1 1=1,15 
TFR3TCn=TSP ' 

Tfn=TEP 

CQI^TIhUE 

ftT»TCl) 

T8=T(A'} 

DO 51 J*l,365 
CALL VI8HA(TES,JJ 
TEW{JJ=TES 
CONTINUE 

PRINT 83, (TE4<( J||J»*.,365,i 
FORMATUX,l6F7.2l 
TEMX=Tgw!lJ 
DO 58 1*2,365 . • 

IT(TEHAX®lT,TEVj{ in TE?«AX*TEWI I ) , 

■CONTINUE 
TF:NIN*T£W(il 
00 61 1*2*365, 

IFCTpiN.GT.TEWUJI Tgl^IN»Te««n 
CONTINUE 
00 63 1=1,365 

IF{TEHIN*lQ*IEWCni TEMD»FLOATCn 
CONTINUE 

A^CTlMAX-nENINJ/a*© 

BB*IT6MAX~T6MINJ/2.':^ 

HS»18*3 
P1*3*14285 
00 49 1=1*365 

TFA*A>“BB*COSU2«0*P1*IFLOAT( IJ-TEHD) 1/365*01 
TeASaJ*TGA 
CONTINUE 
00 2 J»JC,365, 

QTS=HS*CAT-T£ASf JJ) 

OTSS { Jl=QTS 

tCl)*(8*T(2HC#T{lJ“QTS*X3J/Xl 
4T»TCIJ 
.TAVRG=T{1> 




00 3 I«2*N1 

T t U = Ib#T( f-l HB«iTC ni J+C*TC 1 1 )/X2 
TB*{fi*T{tliHC*TINH/Xl 


T(N5»TB 
Jf K*1 

frit ClHTCi H6»6*4 
4K».i 

TAVRe-TAVRG+TlJKl 
COWTINUfi _ ■ _ 


I i I IIS6 /fl8aI?4k 



flfASC4J,@fS 


:TA(Il,l*i 



:lfK 


«'QOC!l} 

Till 


GO TO 99 


99 


I:';: 

ICO 

7? 


J, GO TO 98 
,01 GO 10 77 

CT! I J, I~l,f'0,TEAS{ JItOTS 


46 


57 

10.1 

41 


C.Of\’T!TUl':. 

!F{OlG#iT, 

5 1 = ... 

TFRST(n= 

CCNTI^'UE 
iFCCCOE^PC 
I FtC00fc»Gi 
CQDS=10,O 

loX,4, 

4C=2 
no TC 9a 

rORM/STCiH ,7F9«3J 
FnF.HjST(7F9.3i 
PRINT 4l,I.TR 

PRINT 57,TeW4XtTEMIo,TEMO 

TC€TS1-QTSS( 1>+QTSS<365I 

TaGTS2=0,0 

DO 42 1=2,364,2 

^aC3T S 2=T0CTS 2+QT SS C i J 

CfiCTlNUf; 

T0GTS3«0.a 
00 46 1=3,363,2 
TneTS3=T0«TS3+QTSSC n 
CONTFIUE 

T0CTS=C«T0eTSl+4*^T0QTS2+2,*T0CTS3J#DELTl«i/3» 

’•r;qTS=T001S/365«P 

PUF’-'T 44, TOOTS 

FeRMATClX,*NET HEAT LOSS FROM im SURFACE IN A YEAR=*,Fi 5* 51 
F0RMATax,=fT!:MAX=t,F?.2,4X,*TEMi.'‘-i=*,F7*2,4X,=S‘TEHD=’^,F7,2I . 
F0R-MAT{1K,I3,17F7«2I 


r UKHR \ Iht IS fi.tr Im&l 

FaRMATCiX,*ITR»*,I5I 


C 

$! 


25 



y.9 


*|04»?4S39,19.5: 


,2,11*225,11.982,12*781 


,* o-j ,, jiiU* :j i4t., A Ji.* . 

I.47,16*47,j7*53,i|*64,l9*§2,|l,a6,22*n,23*15t2: 


33 


TOP 

EMC 

SFTC^VIBEA 

SUBROUTINE VIBHACTE,!!) 

DIMENSION TS 1511. F.STSt 511 
DIMENSION POLCilJ 

OIMENStON ATA(3651,ARH(365l,AV{3651,ACFt3651,ABETAl36Sl ,AQSE(3651 
DIMENSION AESTA|365J 

OIMENSION TP(511,flTSPI5l),tRC511,QTSRC5iJ 
CCMMCN ATA,ARH,AV,ACF,ABETA,AqSe,AESTA 
TA»ATA(ni 
RH^ARHlIIl 
^ V=AV{l|) 

CF«ACFfIIl , ' ■ ■ 

BETA»ABETAiIIl 

OSE*AOSiilll . 

FSTA»AE$fACIIl 

Tscii«ia* 

00. S? 1=2,5 
ES 

2* 33, 31*8|.,33»69#,35*65t37* i§, 

■5*3i,§§*33i6|«49,i4-8i3,6|*25,71*87,75«64|79*59f8|*7||S|*i 
4* 54, 102 *10 , Ic*7* 21,112*53, 118*06, 123* 63, 129*85, 136*11, 1^ 

00 33 1=1,51 

FSliasrSf S ( 1 1 

■ QBR«4*14l-eS«{ {273.4Tli4#4-8eTA#C 27i.+TAl«4> 

Q8=19*4*IS*4a40*Q75S*Vl»l‘l£STl-RH*6$TAl 
AOT|«QaR+0S+QH+OE 

I?||]!I1a6TS 

• Si 10|^|»1,51 

iifi 

‘IM«X 
If: 




3*621 
, 26* ’ 
t49./ 


1,52*44,5 


. 0*0 1 so to im 



«T‘SI1I I l-tf ^ lilft; ' 



c 

c 
c 
c 
c 

C“ 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 




•^■cT^RI^Ii^ATICK 
COCLING POND 


OF TW'; i'SATICftL n/"CiRCULATICN li’l THF 


•ois.v 

”01 K I' 


CTION 

CTION 


\Sl IS NOM-DIMENSIQilAL STREAM FU»4CTIu?^ 

ASID IS DIMENSICt^AL STRFAH FUNCTI.Of*^ 

VXN IS ?'.G.‘>DIHENSIO-AL VELOCITY I^: X 
7ZN IS YON-OIMENSIQwAL VELOCITY I?' 2 
STKPX IS STEP SIZE IN X-DIP ECTiO’:. 

STSPZ IS STEP SIZE IN Z-*DIR'XTION 
'•‘X IS m OF OIVISIO':> IN X-DIRFCTI-':^ 

HZ IS NO OF DIVISIO.v IN Z^OIRFCirnw 

IS TOTAL NO OF 6RI0 POINTS I.Y X-ruRCtl I O.J 
7N IS TOTAL NO OF GRID POINTS iM 2“DiPt-XTI0K 
u IS velocity at inllt and out Ls-T 


c 

C INVISCiD FLOy SOLUTION 

C 

SIBJOB 
$IBFTC MAIN 

OINENSIQN ASII(27,21I 
OIMENSION ASIi27*273,ASIOL0C27|271 
' DIMENSION VXN{27*27J,VZN(27, 27!, 

PUAD 50fU,AL,AH,STEPX,STEPZt«X#MZ 
, PRINT 5LCI*AL»AH#STEPX,STCPZ „ :, 

51 - F 0 RHATUHQ,«U=’I',F 8 , ^t^AL«’»‘,F* 3 . 3 »*AH 5 =*tF 8 « 3 ,*STEPX«#,F 6 * 3 »*STlPZ*«» . 


5 

li 


rnir«j p it y * «n I cr A f I 

F0RHATaHQ,«U=’l',F8,3,«L»*,i 
IF6#3J 

C«CAL/AHI**2 
ITRMAX»250e 
00 2 I»2,NZ 
00 Z J*2#MX ■ 

as!ol5<!,41*asici,4I 

•COINTINUg 
NM*MX +1 
NN-NZFl 
DO 3 

i!lliWh2«2o 

CONTINUE 
DO 4 I«2iNN 

CONTlte 
M0=NN~2 
DO 5 I»2,NI 
ASICi,l 1 «l 000«0 
CONT 
lTR=*f 
ITR*’ 

00 6 
00 4 



TR+1 
4*2# NX 
1 * 2 , NZ 


AS 

ICI . , 

CONTINUE 
IG»I3 “ 


1 : 



■?il,41+ASI(I~l,4J+C^ASI(l,4+l}+C#ASIM,4«llJ/C2*Fa** 


fitl#4ll 

F ...,., 








65 


Si-'.MR 

C 

c 

c 

c 

SISUQ 

$iaFT 


on 65 i*l,NN 
Of] 65 

■^SJI(I,J}"ASHI,J)/i>ai5«0 
. Cllf’’T!NUe 

PRINT 14 . , 

FORMATClHCft’^THE VALUE OF OaOCiNrU-lSIONAL STfi£:AM FUNCTION=*',/} 

TO 18 I^l.Nt*! 

PRr-lT 12 ,!aSIICI,JI,J= 1,24J 
CfINTIMUe 

F0RM4T{1HC',24F5-2I 

fnRMAT(5F9.3,2I4) 

STOP 

] m 

¥. ■ ■■ . 

CRggPIWG FLOy SOLUTION 


Oi^EKSlON 
niMEK-SIQH 
nil^ENSlON 
P'.;j 50 sOfU 
PP-INT 51, 


ASIH25,ir^5l 

ASn25,iaf>},ASIOLDC25 


A 3 I I z: 3 , tur^i , AOlULUliiSI, 4 

V X N C 2 5 * 1 0 5 J , V Z N I 2 5 * i 0 f? J 


'i An% £->* 1 9 '?«.!'( \ 

,AL,AH,STEPX,STEPZ, 

U,AL,AH,STHPX,STfiPZ 


:ji&,u,«t.,Ari, 3 icKA» 3 jnr£. 

FORWATClHD,»U=*,Fa, 5 ,*AL«*,F 9 
iF6.3) 

Ca(AL/AH)##2 
CC>=( AL/AHJ*^4 
DA=STSPX’«'«4 
r>B*{ST£PX*STEPZJ’i'#2 
ITPMAX*I400 
OC*STFPZ«*4 
RB»4«i'DC+4**C*D8 
»C*4.«‘0B=^C+4*#CC*DA 
6 £)* 6 .*DC+ 8 «^C*C 8 + 6 ,^CC*DA 

MNfNZ+1 

M*NZ-1 
■TO 2 1*3, N 
DO 2 Ja3,N 

ASI 0 LDCI,j 1 »ASiCI,JI 
CONTINUE 
00 3 J*|,NM 
AS|C1,4|*0*0 
ASI{2i4l*G.Q 
ASIiNZ,4I*S»0 
CONTINUE 
00 4 1*3, N 
AS|CI,«XI*C,0 
CONTINUE 
00 5 1*3, N 
ASI(I,2I«1000*0 

iTRalTR+l 
00 60 1*3, N , 
AS|II,HNl*AS|a,N) 

CONTINUE 
■DO 6i 4*3, N 
■■ ASIiSN,4t-ASI{H,4> 

CONTINUE 
00 6 4*3# N 
DO 6 1*3, N 

!§** I § i f I ■*’ I » 4 1 ■* i 1 1 1 i "" 1 • *11 

, A 6 «Aii|lt 4 '*iH«II!lf 4 “Il 


C25|IC5J 

HX,NZ 

9 . 3 ,*AH**,F&* 3 ,’»STEPX=»,F 6 » 3 ,: 


'STEPZ**, 


"" 1,4 

Tfii 


•IHASlf I-1,4+1HASII1- 


mmm 

■U'AFVmI 


rflilAf-l|jA^"il*C*l 38 «AC* 8 C*AO*CC*OA^If l/SO'' ’ 'U' 


, :.v I s-rf w>,-r 











13 


II 


S3 

21 

22 

21 

24 


25 

50 

12 


AT f?Gr;..CF 


!4J 


$EWTRY 


ccr.ni^ju..: ■ 

IFCBIG.L*’--, 1«2I t;r i i . 

IFC nR«GF»iT?HAXJ 6/ TH I 
DO 7 J«3,^ 

ASIOLD?trJJ = ASICI,lJ 

CONTINUb 
10 TO II 
PRi'n |3»!T8 
F-aRMATllH f^THE NO GF 
no 65 1*1, NN 
DC 65 

■’'SII {I,4)*4SI{I, JJ/lDbr.-F 
CONTINUE 
■^RIMT 14 

FORMAT CIHC),5nHF VALiiF CF qQ"» C I «•■ r M ■.■>*'■,,% L Glh A-’^' FUfXO'I ON* , /I 

na 18 

PRINT i2,UsntI,J)f J*2,2ll 

CONTIsWF ■ 

00 if J='^,FX 

yxNi2,ji»Usnc3,j}"ASiii p/ 

VZNC2,J>=:?."' 

, ■ 

VZNCN2, JSaa*G 
CONTINUE 
DO 32 I=at^Z 
\f im i,2i»o.o .. , 

¥XN( I,MXI«C3.u 
COMTINUS 
00 63 l^aiflZ 
YXWCItZJaCi.Q . 

CONTIMUe . . . 

00 20 Ja3,« 

00 2^ 1=3, N 

VX!SCI,'JJ»CASII( I+l,JI“ASII(I“i,Jn /( 2«*.SfrPZ) : 

C0NTINU2 

DO 11 3=|*Jj 

I,j|2|ASIftl,J-l)-ASIl(l,J+in /Cl.^GTFPX) 

CONTIWE • : 

PR INI 22 '*■ 

FORMAT! VALUE OF K- CONPDNgNT OF VELOCITY*, /I, , 

00 23 1=2 |NZ 

PRINT 12, {VXN!I, 41, 4=2,211 

CONTINUE 

PRINT 24 

FORMAT! |HO,*THE VALUE OF Z-CCMPONFNT OF VEICCITY*,/! 

DO 25 I»2tNZ 

PRINT 12,!VZN!I,4S,J*2,2I1 

CONTINUE 

FORMAT! 5F^*3, 2141 
FORMAT!lH0,20F6,31 
STOP 
FPC 



